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Abstract 

Recently, deeper understanding of QCD emerges from the study of the AdS/CFT 
correspondence. New results include the properties of quark-gluon plasma and the 
confinement /deconfinement phase transition, which are both very important for the 
■ scenario of the QCD phase transition in the early universe. In this paper, we study 

some aspects of how the new results may affect the old calculations of the cosmological 
QCD phase transition, which are mainly based on the studies of perturbative QCD, 
{Sj ■ lattice QCD and the MIT bag model. 
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^ ! 1 Introduction 

Phase transitions can produce relics, affect the anisotropics of the universe, or have other 
observable consequences; hence, it is very important in astrophysics. A particularly impor- 
tant phase transition is the QCD confinement/deconfinement phase transition, in which 
the deconfined quark-gluon plasma (QGP) phase transits to the confined hadronic phase. 
By assuming this phase transition is first order and that it has nonzero surface tension, 
it suffers chronologically the processes of supercooling, reheating, bubble nucleation, and 
may produce relics such as quark nuggets. For the up to date reviews of the cosmological 
QCD phase transition, see [H12]- 

Recently, deeper understanding of QCD emerges from the study of the anti-de Sit- 
ter / confor mal field theory (AdS/CFT) correspondence [3]. In its prototype version, type 
IIB superstring theory on A0IS5 x S 5 is dual to M = 4 U(N C ) super Yang-Mills (SYM) 
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theory in (3 + l)-dimensional spacetime [I]. Generally speaking, conventional quantum 
field theories make sense only in the perturbative regions, where the 't Hooft coupling 
A4 = gYM^c is small; however, the dual gravity theory is easy to handle when the su- 
pergravity (SUGRa) description becomes reliable, that is, in the strong coupling region. 
Hence, we can use the AdS / CFT correspondence to study field theory in the region where 
perturbative approaches are not applicable. It is also believed that the generalization of 
the AdS/CFT correspondence can realize some more realistic systems, such as QCD-like 
theories or QCD itself, which have running coupling constants (hence are not conformal), 
fundamental matter, and reproduce some phase transitions. Recent reviews on the con- 
nection between string theory and QCD can be found in [5j E]. 

New observational results from the relativistic heavy ion collider (RHIC) \J\ data tell us 
that the shear viscosity of the hot plasma is very small [8]; thus the QGP at temperature 
T > T dec should in fact be strongly coupled [HI [10], rather than asymptotic free as we used 
to think about it, where T^ ec is the critical temperature of the confinement /deconfinement 
phase transition. Hence, all phenomenological applications of QGP, which are based on 
perturbative QCD or the MIT bag model [TH [12], should be reconsidered. These appli- 
cations include the neutron stars/quark stars and the cosmological QCD phase transition. 
AdS/CFT provides an excellent tool to study them. Because of its strong interactive 
nature, it can explore the properties of QGP and the confinement /deconfinement phase 
transition, in both the high temperature and high baryon number density regions. How- 
ever, in this paper, we will limit our focus on the property of high temperature region, 
which is important for the cosmological QCD phase transition. 

The organization of this paper is as follows. In Sec. [2, we present the results of QGP 
and the confinement /deconfinement phase transition from AdS/CFT. In Sec. [3j we study 
how these new results affect the conventional scenarios of the cosmological QCD phase 
transition, including the nucleation rate, the supercooling scale and the mean nucleation 
distance. We summarize our results in Sec. [H We will set h = c = k = 1 throughout this 
paper. 



2 The Thermo dynamical and Hydrodynamical Quan- 
tities of QGP Reconsidered 

2.1 Entropy, Free Energy, Energy and Pressure 
2.1.1 Af = 4 SYM theory 

The gauge fields of large- N c Af = 4 SYM theory are described as open strings ending on N c 
Dirichlet 3-branes (D3-branes). In the large 't Hooft coupling limit, A4 3> 1, the entropy 
density can be calculated from the Bekenstein-Hawking entropy [13] of non-extremal D3- 
branes with Ramond-Ramond charge (RR-charge) N c , which is [14J 

s = y Y C 2 T 3 , (1) 
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where T is identified with the Hawking temperature of the black brane. The result is 
only 3/4 to that of the free gas case s = (2tt 2 /3)N 2 T 3 . It was argued that the entropy 
density can also be calculated from the action / by V f = TI = Ve — TVs |15j . thus / = 
— {it 2, /%)N 2 T [IE]. The sound mode dispersion relation of hydrodynamical calculations in 
the strongly coupled limit gives c 2 = dP/de = 1/3 and P + e = Ts [UJ ; hence, both the 
energy density e and the pressure P in the strong coupling case, should be only 3/4 to 
the value of weakly coupled case, which is consistent with the free energy result from the 
action /. In fact, all CFTs' have similar equation of states (EoS's) up to some numerical 
factors [IB], and what we presented above is just a trivial example. 

For the case with not-so-strong coupling, the leading correction is calculated from the 
action J, which reads IT6l 



s 



jj + fc(3)(2A 4 )- 3 / 2 



(2) 



comparing to the weakly coupled case [TjJl [2U] 



So 



2tt : 



:A 4 + 



3 + y/2 



71- 



Af + 



(3) 



The 3/4 factor reveals the intrinsic difference between strongly and weakly coupled system. 



2.1.2 The QCD-like Theories 

However, CFTs are very different from QCD in many aspects. For example, (i) their cou- 
pling constant A4 does not run, hence they experience no conventional phase transitions, 
and (ii) they can only describes fields in the adjoint (color) but not in fundamental (flavor) 
representation of the gauge group. The confinement /deconfinement phase transition is al- 
ways understood as a Hawking-Page phase transition [21] between two background metrics 
with different free energy density / [15] (except the scenario of [22l 1231 [24] ) . The free energy 
density of the system can be calculated from the volume of spacetime J d D Xy/g, and the 
stable spacetime configuration has the lowest /. Flavors are often added by iVf spacetime 
filling (flavor) branes [251 12S] ; however, calculations can be done only in the probe limit 
(exact quenched approximation), Nf <C iV c . Many efforts have been spent to construct a 
more QCD-like dual theory. As a phenomenological discussion of their applications to cos- 
mology here in this paper, we do not want to compare their similarity and dissimilarity in 
detail; however, to make our results more concrete, we do not limit our discussion to some 
special model. We will reveal the bottom-up way (the AdS-QCD approaches) including 
the hard-wall [271 12B] and soft-wall [29] models, the top-down way including the D3-D7 
system [30l EH [321 [33] and the D4-D8-D8 system (the Sakai-Sugimoto model) [3U[35], and 
also some other phenomenological approaches. The comparative theories include the MIT 
bag model [HKI2], the fuzzy bag model [SB] and some lattice results. Most gravity dual 
theories are limited to the large-iV c limit; however, our QCD has N c = 3, which makes 
quantitative applications of the AdS/CFT results difficult. We will try to compare the 
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disagreement between N c — > oo and iV c = 3 by some lattice results [37J. Because of the 
context of this study, we will always assume that the chemical potential \x = in this 
paper, hence the relation between the free energy density and the pressure is / = —p. 

Let us first discuss the AdS/QCD approaches. In the hard-wall model, a cutoff is set 
in the infrared (IR) region to form a slice of AdS^, which makes the boundary theory 
confining J27J 12H]- The two solutions of the Einstein equation are a cutoff thermal AdS 
and a cutoff AdS with a black hole. For the Ricci flat horizon case 1381 [391 
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where the subscript h indicates the confining phase, q indicates the deconfining phase, 
n\ = 8TrGs describes the gravitational coupling scale, and L is the radius of the AdS space. 
For the spherical horizon case with sufficient small IR cutoff tq, we have [10] 
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(5) 



where Q 3 = 2vr 2 and r+ = (3/87rT dcc )( v /9T 2 /T| ec - 8 + 3T/T dcc ). The latter case has little 
physical applications; however, it has thermodynamical properties similar to the soft-wall 
case. 

In the soft-wall model, the IR cutoff is replaced by a smooth cap off, which is realized 
by the dilaton term in the Einstein action |29j. The difference of the free energy density 
of the two phases is [38] 



fq ~ fh 



4r3 



7T 4 L 



rp4 



e~ x (x - 1) + - + x 2 Ei(-x) 



(6) 



where x = (T dcc /0.4917287rT) 2 , and Ei(-x) = - ' e"' /tdt. 

For a ten-dimensional 'AdS/QCD cousin" model with the metric of a deformed AdS§ 
black hole crossing some 5-dimensional compact space [41] , the free energy density is 
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which is related to a entropy density s = sT 3 exp (— T d 2 ec /2T 2 ). This model may be appli- 
cable to QCD for 1.2T dcc < T < 3T dec . It has a good asymptotic behavior limx'^ 00 s oc T 3 
as a four- dimensional thermal system, because the contributions of the Kaluza-Klein 
modes are not taken into account. When T dec <C T, the result coincides the fuzzy bag 



model [36] in pure glue case, which restricts -Bf uzzy 
fq-fh = -{s/^)T\l-TlJT 2 ). 
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There are also some other models, like the one defined by some complex metric in [12] , 
the one include a nontrivial dilaton flow deformation [33], or the MIT bag model itself. 
They all have f q — / n oc T^ cc — T 4 , hence are identical with each other up to an overall 
constant. And in fact, for small supercooling, they are much similar to what in Eq. (Tj0). 

Next, we will discuss the top-down scenarios. In the D3-D7 system [301 EH E21 [33], N c 
coincident D3-branes form an extremal black brane with near horizon geometry AdS§ x 
S 5 , while Nf coincident probe D7-branes fill AdS§ (hence, they also extend along the 
radial direction) and wrap some S 3 inside S 5 . When the D7-branes are separated from 
the D3-branes in S 5 , the chiral symmetry and conformal invariance are broken. When 
the temperature is low, the separation is large enough that the brane tension can avoid 
the D7-branes falling into the black brane, hence the branes are "Minkowski" embedded 
outside the horizon. However, when the temperature is high enough, the gravitational 
attraction of the black brane renders the D7-branes a "black hole" embedding [HI H5] . 
The critical temperature is Tf unc j, where the mesons melt. The multi-valued nature of 
the free energy density makes the phase transition first order. Nevertheless, for massive 
fundamental quarks, it is not the temperature of the confinement/deconfinement phase 
transition, which occurs at some Td ec < Tfund- There is as yet a lack of suitable models of 
confinement/deconfinement phase transition within D3-D7 system. The explicit solutions 
of f(T), s(T) and c s (T) are shown numerically in [HJ [45] . For our purpose, we will not 
discuss this "melting" transition in detail; notwithstanding, we take notice of some of its 
critical parameters which can be compared to that in the confinement/deconfinement phase 
transition. The discontinuity of the entropy density in the phase transition point is 

A /V /V T 3 
As(T = T fund ) ~ 0.066 x -L^It^ ^ 0.032 Km s fund , (8) 

which is proportional to N c Nf, because only the contribution of the fundamental mat- 
ter is taken into account. The entropy density of massless quarks is lim^oo Sfund = 
A 4 iV c iVfT 3 /i6, and the entropy density attributed to gluons is as what in Eq. (pQ). The 
superheating and supercooling ranges are (by the system itself rather than by impurities 
or perturbations) 

T ■ T 
A< = 1 — ~ 0.0019 and A> = - 1 ~ 0.0083. (9) 

Tfund -^fund 

The speed of sound also deviates from l/v3 nontrivially when T approaches Tf unc j. How- 
ever, unless in the extreme supercooling case, c s would not be vanishing. 

Does this "melting" transition happens in QCD? This is an intractable question. Even 
if we neglect the influences of the large- N c and the probe simplifications, we will still need 
QGP remaining strongly coupled at Tf un d; because when it is weakly coupled, the melting 
of the mesons should be a crossover. As we will see later, the numerical values of As, A< 
and A> are all much smaller than the typical confinement/deconfinement case; besides, 
we don't really know how to estimate the surface tension (Xf un d of this phase transition. 
In addition, melting of different mesons may be asynchronous in QCD. Of course, if it is 
indeed a phase transition in QCD, it can also affect the evolution of our universe. 
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In the Sakai-Sugimoto model [Ml [35] , when the temperature is low enough, the N c 
coincident D4-branes are compactified on a supersymmetry-breaking spacelike S 1 to make 
the low energy QCD-like theory (3 + 1) -dimensional, while the Nf D8-D8 pairs (with 
D8 and D8-branes coincide respectively) cross the S 1 circle at some characteristic points. 
Gauge bosons are regarded as massless modes of open strings with both ends on D4- 
branes, while fundamental fermions correspond to open strings with one end in some 
D4-brane and another end in some D8 or D8-brane. However, when the temperature is 
high enough [46J, to make a lower free energy, the compactified D4-brane direction is not 
spacelike but in fact timelike. This is the confinement /deconfinement phase transition, 
because the topological change of the spacetime makes the expectation value of a temporal 
Wilson loops change from (W(C)) = to (W(C)) ^ 0. The spontaneous chiral symmetry 
breaking is understood as when the iVf D8-branes and iVf D8-branes merge at some radial 
position uq away from the horizon (where we live), which happens at some temperature 
higher or equal to T^ ec . The difference between the free energy densities of the two phases 
can be calculated from the DBI action. This phase transition is first order, 

409607r 11 / s (^ s iV c )iV c 2 

Jq-Jh- ™ 7f K 1 - J dcJ- l iU J 

I ZM 1 dec 

For its AdS e non-critical string "cousin" model [47], f q — fh oc — iV c 2 (T 5 — T| ec ). In the 
Sakai-Sugimoto model, one always have the speed of sound c s = l/v 7 ^ [IE]. Because 
l s (g s N c ) = g 2 5 N c /(27T) 2 = g 2 N c / '(2vr) 3 T dcc = A 4 /(27r) 3 T dec , and A 4 iV c /2167r 3 ~ 7.45 x 10~ 3 
from meson spectrum [35] , we see that the coefficient of Eq. (fTOj) is really huge. However, 
these results are quantitatively far from QCD; the unwanted Kaluza-Klein modes of the 
compactified dimensions cause the theories lacking of the asymptotic UV behavior / oc 
N 2 T A while T —>■ oo. As we will scale all these theoretical models to QCD by their high 
temperature behavior, we will not consider the Sakai-Sugimoto model from now on. 

There are also some more phenomenological approaches to the EoS's of the QCD-like 
theories. Giirsoy et al. considered a five- dimensional gravity theory coupled to a dilaton 
field [491 [50]. The thermodynamics of this system can be determined uniquely by a positive 
and monotonic potential V{\) = 12[1 + A + Vi\ 2 ® log p (l + V2A 2 )], where <p = log A is the 
dilaton field [SU [52] . The theory is confined when Q = 2/3 and P > 0, oi Q > 2/3. 
After chosen some specific potential, the temperature is fixed uniquely by the horizon 
value of A, and the EoS can be given by some numerical calculations of the black hole 
configuration while varying A(r#). The aim of this model is still limited to explain the 
finite temperature large- N c Yang-Mills theory by these authors; however, we may expect 
that can tell us something more about QCD. 

Gubser et al. considered another five- dimensional gravity theory coupled to a single 
scalar [22]. Based on a lot of assumptions, it is shown that the potential of a scalar field 
V((f)) and the EoS of the boundary theory have one-to-one correspondence. The results 
may be applicable for regions both > and < T dec . Various V(0)'s correspond to different 
EoS's, include crossover, first order, and second order phase transitions; hence, the authors 
expect their model can mimic the EoS of QCD. We choose in this study the potential 
V(cj)) = [-12 cosh (70 + b(f) 2 )]/L 2 with 7 = a/7/12 and b = 2 for the first order case when 
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making comparison with other models □, but keep 7 as a free parameter in the discussion 
of Sec. m 

The phenomenological model in [49l [50] coupled to the dilaton potential V(cf>) has a 
more solid theoretical foundation; however, the calculation of the EoS's is more complicated 
than the latter one. As we need to exploit a whole family of EoS's for our astrophysical 
purpose (especially in Sec. 13.3.21) . we will limit our discussion to the latter model. It should 
be noteworthy to review the astrophysical application of the first model, especially after 
some quantitative comparisons between it and the lattice results that have been done. 

The comparisons of the free energy density /, the entropy density s and the square 
of sound speed for various models, are shown in Fig. CD, E] and [31 We scale all ther- 
modynamical quantities by T dec and lim T ^ 00 (»)/(») qi SB = 3/4, where "SB" denotes the 
Stefan-Boltzmann values of thermal quantities in the corresponding QGP phase, except 
the model discussed in [511 [52]. The scaling relation is based on Eq. ([1]) and the fact 
that all gravity dual theories are strongly coupled; we assume that all field considered are 
UV conformal, and the coefficient 3/4 is universally applicable for them all. The model 
in [511 E2] is excluded, because it is indeed weakly coupled in the UV region and asymptot- 
ically Stefan-Boltzmann. The rescaling is of course reasonable for the entropy s in Fig. [21 
because in the large- N c theories, s q oc N% and oc AT?; hence the latter one can always 
be neglected. For the free energy /, things are a little more subtle, since the UV cutoff 
introduced by the computation [15] do not ensure /h oc N®. We make the statement as an 
assumption by using some appropriate counter-terms. Although some models (e.g., [UJ 
and [22]) aim directly at QCD itself, we assume the superheating contributions of /h and 
in their models can also be neglected comparing to QGP while T — > 00. The numbers for 
the classified models are tagged in Fig. [lj In Tab. [IJ we list the maximal superheating and 
supercooling scale A> and A< for various confinement /deconfinement models, and also 
the "melting" transition in [HI [45]. The existence of A> and A< indicates a completely 
different phase transition process comparing to the old one; for the range of superheat- 
ing or supercooling is no longer caused by impurities or perturbations, but caused by the 
theoretical system itself. It can be seen that the "melting" values are much smaller than 
the confinement/deconfinement case. In Fig. HI we compare the latent heat from the 
theoretical models list above, and from the lattice calculations for various N c > 3. 

Heuristically, we see that for temperature 2T dec < T < 4T dec , models (2), (3) and (4), 
(5), (7) look similar to each other; and for temperature T ~ T dec+ , models (2), (3), (6), (7) 
look similar. Models (2), (3), (6), and (7) have the latent heat a little too small compared to 
the lattice result of the large- N c theories, and the latent heat of model (1) and the original 
MIT bag model seem too large. The discussions of models (6) and (7) may be a little more 
unreasonable, because the free parameters in the dilaton potential are chosen arbitrarily. 
The divergence of these models rise because they are all quantitatively far from (large- N c ) 
QCD. To avoid the unnecessarily complicated details of these models in our discussions in 

1 Our EoS's agree with [22] qualitatively, and can catch all the expected limits; however, it has some 
small quantitative divergences rising from the numerical trickings. These divergences can affect the critical 
parameters we choose for the second order phase transition, but will not affect the main results of our 
discussions. 
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Table 1: Comparison of superheating scale A> = T max /Td ec — 1 and supercooling scale 
A< = 1 — T min /Td oc for various models. The models are numbered as in Fig. [IJ The A> 
and A< in the "melting" transition [HII15] are simply replaced T dec by T fund . 



Model No. 


A> 


A< 


1 


oo 


0.159 (oo) 


2 


oo 


0.057 


3 


oo 


1. 


4 


oo 


1. 


5 


oo 


1. 


6 


oo 


0.111 


7 


0.046 


0.011 


"melting" (T fund ) 


0.0083 


0.0019 



Sec. El it is worthwhile to ask what kind of feature they have in common. We argue that 
(i) the EoS of real QCD should be softer than the bag model, and (ii) there exist some 
intrinsic maximum supercooling scale to be achieved, in contrary to the old belief that the 
range of supercooling is caused by impurities or perturbations. 

It is interesting to argue in what conditions the bag-model-like theories can still be 
applicable. The renormalization-group-improved perturbation expansion method tells us 
that, when the strong coupling constant a s increases, the bag constant -Bmit decreases [53] . 
Although this result is only suitable for the perturbative and zero temperature regions, it 
suggests us to treat the bag model carefully. However, there are indeed a lot of gravity dual 
theories whose boundary field theories have bag-model-like thermodynamics [381 139| H2 | |4"3] . 
which do well for explaining meson spectrum or other physical applications. 



2.2 Shear Viscosity and Bulk Viscosity 
2.2.1 M = 4 SYM theory 

The shear viscosity of large- N c M = 4 SYM theory in the large 't Hooft coupling limit, 
can be calculated via the Kubo relations. The result is [551 \57\ 
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And for bulk viscosity, conformal property requires £ = 0. It was argued that this value is 
always available for theories with holographically dual supergravity descriptions 



For the case of large but finite 't Hooft coupling A4, we have [60, 61J £ 

-3/2 
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which can be compared with the weakly coupled case 

77 6.174 



s Xj In (2.36/ VXD 



(13) 



2.2.2 The QCD-like Theories 



As for the thermodynamical case, people follow the top-down and bottom-up routes to 
discuss the hydrodynamical quantities of QCD-like theories. However, there is a lack 
of lattice results to be compared with, because lattice QCD is incapable for real-time 
behaviors. 

To break the conformal behavior of AdS/CFT, one easy way is to consider Dp-branes. 
The result is [63] c 2 = (5 — p)/(9 — p), 

V 1 A C 2(3 -p) , 1A . 

- = — and - = —. ^. (14) 

s 4i 77 p(9 — p) 

For the case of compactified Dp-brans, the relations for c s and rj/s are the same as before, 
but the relation for (/r] has to be modified to 

C 8d-2(9-j>)(<i-l) n (\ a 



d(9-p, 2 {d- C -)' <15) 

which is consistent with the Sakai-Sugimoto model's result (/rj = 4/15 [48] for p = 4 and 
d = 3. 

To take into account the contributions of fundamental matter, one can consider the 
D3-D7 system. The result is 



77 = -N 2 T 3 



i+ii^V— ) + 

8vr 2 iV c V^q 



(16) 



where M q is the quark mass, h(x) is some smooth function connects h(0) = and h(oo) = 1 
by a crossover around x ~ 1, with the entropy density s = (n 2 /2)N^T 3 + Sf un( j already 
been discussed in Sec. 12.1.21 Similar calculations for the Dp-Dq-Dq system including the 
Sakai-Sugimoto model can also be done. 

For the models of five-dimensional gravity coupled to some dilaton fields, the bulk 
viscosity can be calculated directly by the Kubo formula [231 121] • C can De estimated by 
the numerical solution of the metric. 

Based on the discussions above and also some other evidences, people conjecture that 
there may be some universal bounds of shear viscosity 77/5 > 1/Att (or h/A^k-Q when getting 
back the units; also called the Kovtun-Son-Starinets (KSS) bound) for all physical systems 
in Nature [58| |59 | [65]. and of bulk viscosity C/77 > 2(1/ p — c 2 ) for theories with holographi- 
cally dual supergravity descriptions [66J . The universality of these bounds suggests that we 
can use them as critical parameters for the properties of QGP; however, different opinions of 
them exist in literatures. Clues from the generalization of the second law of thermodynam- 
ics (GSL) suggests some origin of the KSS bound from very basic physical principle 



9 



nevertheless, various theoretical models have being constructed which violate the bound, 
both from quantum field theory [681 ESI EDI HQ and from AdS/CFT itself [721 [73]. Fortu- 
nately, the latter violation only loosens the bound a little, to rj/s > (16/25) (l/47r), for the 
constraint of causality [731 [71]. I n addition, using the model constructed in [231 El] to cal- 
culate the bulk viscosity of the potential V(<p) = [—12 cosh (jcj) + b(f> 2 )]/L 2 , can sometimes 
violate the bound given in [66J. 

For concreteness, we go back to the case of QGP itself. Let us first discuss the shear 
viscosity rj. Although some theoretical arguments suggest us that rj/s should be much 
larger (maybe by a constant of ~ 7) than 1/Att in the strong 't Hooft coupling limit, 
because it is much larger than the M = 4 SYM theory case in the weak coupling limit [62] , 
RHIC results tell us that the rj/s of QGP nearly saturates [HI [751 US], or maybe even 
violates [76] the KSS bound. 

There are few discussions about the dependence of parameter 77 on the temperature T. 
It has been done in the hard-wall and the "AdS/QCD cousin" models [77] ; nevertheless, 
they both always have rj/s < 1/47T, which violate the KSS bound. Naively, one can estimate 
it by some phenomenological relation 

77 ~ elcs, (17) 

where I is the correlation length; however, it is very hard to make quantitative computations 
by this formula. Some interpolation between strong and weak coupling regions may be also 
possible [78], as the perturbative QCD result of 77 in the weak coupling region is rather 
credible [79J. 

For the case of the bulk viscosity (, lattice results of gluodynamics show that it rises 
sharply when T — > (Td ec )+ [EH EH [82] , which are qualitatively consistent with the fact that 
c s drops there. Although ( cannot be calculated in the supercooling region T < T^ cc within 
the lattice framework, we assume from AdS/CFT that it varies smoothly while cross the 
phase transition point. 

2.3 Surface Tension 

Very few works exist addressing the surface tension of the confinement/deconfinement 
phase transition from the AdS/CFT viewpoint. For this purpose, two separate metrics 
with different topologies, both have (3 + l)-dimensional translational invariance within 
"our world" (as assumed by all the models in Sec. I2.1.2p . are not suitable; as we need 
nontrivial metric change along the direction of "our world" . Some relative discussions can 
be found in [S3]. Deconfined regions map to some pancake-like black hole solutions, whose 
interior resembles black brane; however, they have domain-wall-like boundary to smoothly 
connect with the confined gravity solution. Hence, the hadronization of the plasma balls 
can be understood as the Hawking radiation of the dual black holes. Although this work 
aims particularly at the large- N c gauge theories, some other authors believe that dual black 
holes are in fact produced inside of RHIC [SU 135] . 

The concrete calculation is based on some finite temperature Scherk-Schwarz compact- 
ificational metrics, which have covering space asymptotically AdSd+2 near the boundary. 
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Both the time direction r and a spacelike direction 6 are compactified to some circles S 1 ; 
however, the 9 circle shrinks to zero at some finite u = Uq in the confined phase, rather 
than the r circle shrinks to zero in the deconfined phase. The metric of the domain-wall- 
like boundary can be solved numerically, and the surface tension can be estimated by it. 
The surface tension a is rounded to numbers 2.0 e q (Tdec)/^dec for d = 3 (a hence (2 + 1)- 
dimensional gauge theory) and 1.7 e q (T dec )/T dec for d — 4 (a hence (3 + l)-dimensional 
gauge theory), a oc e q oc N% is a natural result of the scaling of the classical gravity action. 
The aftermath of this fact is discussed in Sec J3. 3.1"! 

However, there are some relevant discussions of the surface tension a, based on both 
lattice gauge theory and the MIT bag model. The lattice results of a for the pure gluon 
SU(3) gauge theory are around 0.02 T^ ec [861 E0 EZ]. In the MIT bag model, the contri- 
bution of a is divided to an intrinsic and a dynamical surface tension [55]. The intrinsic 
surface tension 07 is suggested to be very small; however, we do not know how to calculate 
it in this framework. The dynamical surface tension <jd raises from the modification of 
the fermion density in the phase transition surface; hence, it depends sensitively on the 
strange quark mass. Detailed calculation shows that od is at most (60 MeV) 3 [88]. Notice 
that the bag model results are only valid for the zero temperature case, and the lattice 
results do not consider fundamental quarks (which is supposed to be crucial in the bag 
model discussions). However, these results may suggest that a is not very large. 

3 The Cosmological QCD Phase Transition Reconsid- 
ered 

If the QCD confinement /deconfinement phase transition is first order, just as what the 
application of a Hawking-Page phase transition indicates, our universe underwent that 
transition when it was about 10~ 5 s old. Generically, if the surface tension of the transition 
interface is nonzero, the universe should be supercooled for some scale before nucleation 
indeed happens [EH [90]. After the supercooling stage, some hadronic bubbles are created; 
they may then expand rapidly as both the detonation [HU E21 E3] and deflagration [921 [9U 
193] waves. For the deflagration wave case, the latent heat released by the phase transition, 
reheats our universe back to T& cc . After that, the phase transition goes along synchronously 
while the universe expands, and converts the denser QGP matter to the less-dense hadronic 
matter mildly. The mean distance between the hadronic bubbles, is calculated in [951 ISE 
[nZl IHHl IHH] • After about half of the QGP matter has been converted, the hadronic bubbles 
are replaced by the QGP bubbles. As the phase transition goes on, the QGP bubbles 
disappear more and more rapidly [100] . Baryons may be concentrated in the QGP bubble, 
and relics such as quark nuggets may be produced [101] . Some panoramic description of 
this phase transition can be found in |102j . and some up to date review articles are in [H[2]. 

The process we described above is called homogeneous nucleation. We will not con- 
sider other possibilities such as heterogeneous nucleation [961 EH] or inhomogeneous nucle- 
ation |103] in this paper, because they are less sensitive to the intrinsic properties of QCD 
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(hence, less sensitive to the AdS/CFT results) than the homogeneous case. In addition, 
we will not consider the late stage issues of this phase transition, such as the stability of 
quark nuggets, because the zero chemical potential assumption is no longer suitable there. 
We leave the relative discussions in the follow-up studies, by which the results from the 
finite chemical potential AdS/CFT correspondence can be used directly. 

3.1 The Nucleation Rate 

The nucleation rate of the hadronic phase out of the QGP phase can be calculated as 
in [TIM] 

r = f fi e- Ai W /T , (18) 



where 



T 2 {s q - s h ) 2 Ri 



is the dynamical prefactor to describe the dissipation effect, 



2 /a\W(R 



l 



is the statistical prefactor, and 

AF(«.) = ^J^jy (21) 

3 (/q-Jh) 2 

is the additional free energy of a hadronic bubble of the critical size R* = 2cr/(/ q — /h) 
within the QGP phase, £ q is the correlation length in the QGP phase. For the case of zero 
chemical potential, we have / q — /h = Ph — P q and the enthalpy density lj = sT. 

The prefactor (k/27t)Qo in the nucleation rate formula for various models, is shown in 
Fig. [5j The most important step is how to map the various thermodynamical quantities of 
large- N c theories from AdS/CFT models to real QCD. Our strategy is linearly map («) qi sb 
to the corresponding quantities of the g q = 37 + 14.25 ideal gas model, and map the / q = /h 
and s q = Sh horizontal lines in Fig. [TJ and [2] to the g h = 3 + 14.25 ideal gas model, where g q 
and are the degrees of freedom of the real world at T = T dec ~ 192 MeV [105j before and 
after the confinement /deconfinement phase transition. The coefficient 14.25, contributed 
by the leptons and photons, is almost irrelevant to our follow up discussions, beside the 
ones using the Friedmann equations to describe the expanding universe; hence, we will not 
discuss its rationality. However, the contribution 3 from the pions, actually needs to be 
studied more carefully. Pionic freedom is caused by the fundamental quarks, while 21 of 
37 in g q is caused by the fundamental quarks as well. As nearly all our models of EoS's 
are dominated by gluodynamics, and the contribution to the latent heat or the surface 
tension a by gluons and quarks cannot be discussed separately, this manipulation is in fact 
untenable. However, the quenched lattice method faces the same problem. Nevertheless, 
we take the whole EoS's to describe the thermal quantities in different temperatures, rather 
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than some characteristic parameters like or a. For some models with free parameters 
like in [221 1221 [21] (which we will discuss especially in Sec. I3.3.2p . we may expect that 
suitable choice of parameters can absorb the contribution of fundamental quarks. Hence, 
we expect the calculations below can still reveal some aspects of real QCD. 

As seen from Fig. [51 the strongly coupled nature of QGP can lower the prefactor 
[k/2t[)VLq a lot, mainly by the reason that it has a relatively smaller shear viscosity rj q = 
s q /47r. It is artificial that the lattice results seems much larger than what is in all of our 
models tagged by numbers. The reason is that, the value — 1.4 T^ ec is calculated by 
gluodynamics, but it has been shared naively to both gauge and fundamental particles 
by our simple mapping. As the lattice results indicate, the latent heat of QCD with 
physical quarks may be smaller than pure gauge case, the prefactor may be enhanced. 
The increasing of (k/27t)Q for some not- very-small supercooling for our models is very 
interesting. Beside the reason we erase all the reductions for small supercooling, the main 
reason is when the EoS is not bag-model-like, the latent heat is not as large as in T^ ec 
while the supercooling is large. This can be seen roughly from Fig. [2] and the relation 
L h = (4/3)T(s q -s h ). 

3.2 The Supercooling Scale and the Mean Nucleation Distance 

To estimate the supercooling scale quantitatively, we have some separate criteria. If the 
supercooling is required to complete the phase transition, we need at least one nucleating 
bubble per Hubble volume; that is, T > l/d^At for the Hubble radius dn = a/ a = 
a/ 45/47i" 3 M p i • gq 2 T~ 2 and the nucleating duration At. We may relax At to the Hubble 
time c?h/2, because the resulting supercooling scale is in fact insensitive to this parameter. 
Hence, the supercooling scale can be roughly estimated by T ~ 1/^h- 

To estimate the supercooling scale more accurately, let us consider the deflagration bub- 
ble scenario. The applicable parameter space of this scenario is discussed in [93]. Assuming 
that a hadronic bubble created in the supercooling QGP phase expands deflagratingly, a 
shock wave with velocity v s ^ > c s preheats the QGP matter to stop the new nucleating 
processes there, and a deflagration wave with relatively slow velocity t>d e f burns the QGP 
matter to hadronic matter behind it [921 [9U [93] . The velocities v s h and t>d e f are calculated 
accurately in [93]. The weakly and electromagnetically interacting particles can affect 
these velocities |1071 1108"] ; however, deflagration happens only during the early stages for 
the small supercooling case, when their influences are negligible. When most of the space 
has been swept by the shock wave, the supercooling process ceases. The fraction of space 
which has already been swept by the shock wave is calculated foremost in |109^ 1110] . For 
our purpose, we can neglect the expanding of the universe in the supercooling timescale. 
Hence, the criterion of the supercooling scale Tf is roughly [90] 

An r tl 

— Tv^U - t)*dt ~ 1, (22) 
where tdec(^f) is the age of the universe at temperature T^ ec (T{). This integral equation can 
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be solved approximately by 



V^G ■ e,T(gq + ej djAF/T) 
de q / dT dT 



In (^-tt I c;;,r 

\Z7T 



AF/T 



(23) 

T f 



in which we deal with the Friedmann equations without any assumption about the EoS of 
the QGP phase. 

The numerical result of A = 1 — Tf/Td ec depends on various surface tension a for 
various models, is shown in Fig. [61 For small a, the system follows nicely to the relation 
A oc cr 3 / 2 /L h [96] for fixed L h ; however, when a is large enough, these lines tilt up. 
One reason for these departures from a 3 ^ 2 can be seen from the reduction of Eq. 



for some EoS's with constant L n , which gives A oc a 3 / 2 / a/171 — 4 In (/3/cr 3 / 2 ) for some 
explicitly written positive f3 [90]. The other reason is the effective latent heat released 
drops for some not- very-small supercooling scale for the more realistic EoS's. Nevertheless, 
comparing to the tilting up of d nnc seen form Fig. 0, the effects here for A is really weak. 
Lines in that figure cannot be extended to larger A, in where d(AF/T)/dT — > and our 
approximation becomes inapplicable. In addition, A is totally insensitive to the prefactor 
in the right hand side of Eq. (|23"1) . such as the shear viscosity rj q or the shock viscosity t> sn . 

In a more accurate (and also more sophisticated) way, supercooling scale can be calcu- 
lated dynamically from the time evolution of the temperature |lllj . We don't calculate the 
time-dependent solutions here, because our qualitative QCD theories still have too many 
free parameters, thus intrinsic discussions are not very easy. Nevertheless, we think that 
there should be some interesting results in the not-very-small supercooling regions. 

The mean nucleation distance of the hadronic bubbles in the phase transition era, can 
be estimated by d nuc ~ n(tf)" 1 / 3 and the bubble number density calculated in [1091 1110] . 
Some suitable reductions give d nuc ~ (87r) 1 / 3 w sh /(-d(AF/T)/cit| tf ) [28]. Considering some 
special EoS, we have 



"nuc — 



-57T 



iV3 



u S h de q /dT dT 



e y 2 (p q + e q )rf(AF/T) 



(24) 

If 

The numerical result of d nuc is shown in Fig. [71 It can be seen that for small a, d QVLC cx 
cr 3 / 2 /L h for fixed L h , as is estimated in [96, ,99]; however, when a becomes large, d nuc tilts 
up caused by both a more accurate treatment of supercooling and the drop of L h for some 
more realistic EoS's. Although models (2) and (7) both have some maximum a where A< 
is saturated, their behavior are completely different. In model (2), L n — ► hence d nuc — > oo 
while A — > A<; but in model (7), L n ^ hence d nnc is finite. 



3.3 The Supercooling Scale and the Mean Nucleation Distance 
Once More 

It may not be plausible to consider the dependence of the supercooling scale A and the mean 
nucleation distance d nnc on the surface tension a. New phenomena deviating from the rough 
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analytic estimations A oc cr 3 / 2 /Lh and d nuc oc cx 3//2 /Lh [96], always appear in the regions 
where a is large enough. Although a is indeed a free parameter since we do not know its 
value, it should not be very large both from the lattice results of gluodynamics [861 [83, [37] 
and some theoretical estimations based on the MIT bag model [55, 88J. This issue has 
already been discussed in Sec. 12.31 

Notwithstanding, we can still do some qualitative or quantitative estimations, and give 
some constraints on both the surface tension a and the latent heat L^. 

3.3.1 The Global Constraint of the Surface Tension on the Large- iV c Theories 

In [37] , the authors argued one cannot distinguish the scaling of the surface tension a oc N c 
or a oc iV 2 from their lattice analyses of SU(N C ) gauge theories. However, for the reason 
that we definitely know the latent heat oc iV 2 for a first order phase transition, if 
this transition indeed exists, to avoid a zero nucleation rate in Eq. (Tl8|) . we need at most 
aociV c 4/3 . 

If in some large- N c theories, a dependents on N c sharper than N^ 3 , we can equivalently 
give an upper limit for N c . For the finite temperature Scherk-Schwarz compactification 
model, the domain wall tension a oc e q /Td ec oc iV 2 has been calculated numerically [HS] 
for the compactified AdS^ and AdS§ soliton solutions. Hence, given an explicit expanding 
universe, we can restrict N c by the phase transition happened there. A special example 
to constrain N c of the large- N c CFT in the holographic Randall-Sundrum (RS) I model, 
is given in |112[ I113j . despite of the fact that the concept of the surface tension does 
not intervene their discussions. The exponential suppressive factor in the nucleation rate 
formula, is given by the Euclidean action which has a minimum at T = 1/^3 T c for some 
transition happens at T c . The comparison between the holographic RS I phase transition 
and our model based on AdS/CFT, is given in Sec. [H 

3.3.2 The Extremely Weakly First Order Confinement/Deconfinement Phase 
transition? 

The order of the confinement /deconfinement phase transition for QCD with physical quark 
masses, is still being debated. The lattice results of quenched QCD suggest that it is at 
most weakly first order [57] . However, adding massive quarks seems to make the transition 
weaker, or even gradually changing it to a rapid crossover [114^ 1115] . Hence, one possibility 
to be considered is the extremely weakly first order case. We still assume the bubbles 
expand deflagratingly in this case. 

Naively, both the supercooling scale A and the mean nucleation distance d nuc increase 
reciprocally while the latent heat Lh decreases, base on the rough analytic estimations 
A oc cr 3 / 2 /Lh and d nuc oc cr 3 / 2 /Lh [96]. However, more abundant phenomena can happen 
for more realistic EoS's of QCD. 

These phenomena are caused mainly by two reasons, (i) If the EoS's possesses the 
weakly first order phase transitions, the effective L h decreases when the supercooling scale 
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becomes large. This can easily be seen from Fig. [2] and the relation = (4/3)T(s q — Sh). 
(ii) As a universal property of the Hawking-Page phase transition [21] , there is a minimum 
temperature T m i n < Td ec below which the high temperature phase cannot exist. It is 
illustrated in Fig. El The qualitative effect of the first reason has already been discussed 
in [HE]- We will give here both quantitative effects of (i) for some specific EoS's, and also 
some qualitative effects of (ii). 

For our discussions, we will use the mimicking model of Gubser et al. [221 [23], [23]. 
The reason is that, it is convenient to use its potential V(4>) to construct a first order 
phase transition with decreasing Lh, which then transforms smoothly to a rapid crossover. 
Another phenomenological model including a dilaton field given in [5H [52] may also be 
used, as it has a more solid theoretical foundation. We omit the discussions of it here, 
because the work for this model itself is still on its way, and the calculation of the EoS's is 
more complicated than the former one. Some qualitative properties, such as A< decreases 
with decreasing L h , are supposed to be universal. 

As the potential of the Gubser et al. model V(<f>) = [—12 cosh (70) + b<p 2 ]/L 2 has two 
parameters 7 and b, in fact, our method applies to a wide range of models (potentials) 
with one free parameter. We fix b = 2 and evaluate 7 G [0.722, 0.790] (formerly we used 
7 = -y/7/12 ~ 0.764); when doing this, the latent heat Lh varies from 0.69 to 3.77 T^ ec . 
The dependence of d nuc and the supercooling scale A on L h are shown in Fig. [9] and [KB 
Observing from the figures, when Lh is large enough, it follows the scaling law A oc L^ 1 
and d nuc oc L^ -1 ; however, for smaller L h , d nuc tilts up because the effective L h drops for 
reason (i). In a large acceptable parameter space, d nnc is not as small as people used to 
think as about ~ 2 cm [99] for the homogeneous nucleation case. For some definite a, 



there exist some minimum Lh,<, where the maximum supercooling A< is achieved. 

What happens if the realistic Lh is smaller than Lh ) <(cr)? Maybe this situation never 
happens in a consistent world. In despite of that, as a lack of the complete origin of the 
surface tension, we just treat Lh and a as free parameters. If this happens, we have the 
bubble number density 



n(t < 



K 



-a 



-AF/T 



2tt V -d(AF/T)/dt 



-d{AF/T)/dtf 



t=t< 



sh 



(25) 



t=t < 



comparing with Eq. f[24|) and the discussions in 



where t< is the time when the mini- 



mum temperature (1 — A<)Td ec is achieved. Because of the exponential suppressed factor 
exp (—AF/T), this situation will lead to a much smaller bubble number density n hence 
a much larger d nnc . One may think that the larger d nuc can help surviving the quark 
nuggets, or provide the inhomogeneous initial conditions of the big-bang nucleosynthesis 
(BBN). However, this scenario is in fact rather hard to appear. We also show in Fig. [10] 
the criterion T ~ I/^hj that is, the supercooling scale needed for d nuc ~ g?h- Because d nuc 
varies too sensitively to the supercooling scale, the corresponding L h has some value very 
close to Lh,<. Hence, to get an appropriate d nuc for our universe, we need to fine-tune L h 
in a very small region a little smaller than <, which is unlikely to be so. 
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4 Discussion and Conclusion 



In this paper, we discussed some implication of the new AdS/CFT results to the cosmo- 
logical QCD confinement /deconfmement phase transition. We limit our discussion to the 
homogeneous nucleation case. The values of the hydro dynamical quantities, like the shear 
viscosity r\ or the bulk viscosity (, can significantly lower the prefactor (k/2it)Q of the nu- 
cleation rate formula compared to the old estimations; however, they can hardly affect other 
characteristic parameters of this process, such as the supercooling scale A = 1 — Tf/Td ec 
or the main nucleation distance d nuc . The new EoS's, which differ from the MIT bag 
model, can affect the phase transition scenario mainly in two ways, (i) As most of these 
EoS's are comparatively more weakly first order than the bag model, it is not adequate 
to treat their latent heat L\ l as a constant. For some not- very-small supercooling, the 
effective latent heat is always much smaller. Hence, d nnc enhances comparing to the old 
estimation d nuc oc a 3 ^ 2 /Lh [96] when a becomes larger or becomes smaller. In a large 
acceptable parameter space of a and Lh, c^nuc is not as small as people used to think as 
about ~ 2 cm [99] for the homogeneous nucleation case, (ii) The high temperature phase 
should have an intrinsic maximum supercooling scale A< based on a Hawking-Page type 
phase transition. This is in contrast with the old belief that the range of supercooling is 
caused by impurities or perturbations. We discussed the possibility that this maximum su- 
percooling scale is saturated in the cosmological QCD phase transition, which may happen 
when this phase transition is extremely weakly first order. If it happens, the nucleation 
distance d nuc can be increased tremendously. However, it is unlikely to be so; because to 
get an appropriate d nvLC for our universe (that is, to help understand the surviving of the 
quark nuggets, or to get the appropriate initial conditions of the BBN), needs to be 
fine tuned. 

Some related works are listed as below for comparison reasons. The nucleation rate 
and also some of its cosmological applications, base on the holographic RS I model, are 
discussed in |112l 1113] . In this model, a "Planck brane" and a "TeV brane" are added to 
the AdSz, x S 5 spacetime with a dual CFT. The "Planck brane" makes a UV cutoff hence 
adds a (3 + l)-dimensional gravity; the "TeV brane" makes an IR cutoff, and the standard 
model fields in it are understood as bound states out of the strong interacting CFT |116j . 
When at finite temperature, to make a lower free energy, the low temperature phase is as 
in the RS I model, but the high temperature phase favors an AdS-Schwarzschild solution 
(duals to the free CFT gas); hence, our universe should suffer a phase transition at some 
T c lower than the Fermi scale. To ensure that the phase transition is completed thus for 
avoiding an empty universe, we need a strong upper bound for N c of the dual CFT field. 
This model has already been discussed in Sec. 13.3. 1[ where we pointed out that an upper 
limit of N c may be universal for some large- N c theories which suffer some phase transitions. 

The phase transition of an AdS/CFT model, in which a (2 + l)-dimensional field theory 
is dual to some (confined) AdS soliton or some (deconfined) black 3-brane metric compact- 
ified in a brane dimension, is discussed in |117] . The supercooling and the rapid reheating 
(hadronization) after it, are considered. Notwithstanding, in the large- N c limit, the slowly 
hadronized phase at the temperature Td ec do not happen in their model. To begin at some 
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supercooling temperature T > 0, the residual deconfined regions after the phase transition 
still hold the energy portion larger than 1/4. In that model, the supercooling scale is given 
by hand, and a lower limit To = (A = 1) is considered. Comparing to that work, what 
we do in this paper is calculating A explicitly within some physical environments (what 
we use is the cosmological QCD phase transition). We use some AdS/CFT models more 
pertinent to the (3 + l)-dimensional QCD than theirs. 

In addition, an interesting relation between the KSS bound and strange quark stars, 
is shown in [118J. The authors argued that, the surface of quark stars at the temperature 
T ~ 80 MeV, has already saturated the KSS bound. 

The question which parallels to the topic we discussed in this paper, is how the RHIC 
results of strong interacting QGP and the AdS/CFT correspondence can affect the research 
of neutron stars and quark stars. The difference is that the deconfined QGP in quark stars is 
mainly caused by its high chemical potential, rather than caused by their high temperature 
in RHIC or the early universe. A lot of AdS/CFT models for finite chemical potential have 
already been constructed; although just as in the finite temperature case, they are mainly 
studied in the large- N c limit. We will leave these issues to the follow-up studies. 
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Figure 1: The free energy density for various models compares to the free gas case. For 
the reason that all gravity dual theories are strongly coupled, we scale (•)/(«) q , sb to 3/4 
as T — > oo, based on the M = 4 SYM result in Eq. (JI|), where (•) can be replaced by any 
thermodynamical quantities, such as entropy, free energy, energy or pressure. The MIT 
bag model [HI [12] and the fuzzy bag model [36] are also scaled to 3/4 by some comparison 
reasons; they can be easily transform back to their original form if needed. For clarity, 
we classify and tag our models by numbers. From the arrow direction marked in this 
figure, the thick lines are for the models (1— >4^5^7^2^3^6) respectively. 
Line(Model) (1) denotes the hard-wall model with the Ricci flat horizon calculated in 
Eq. (T4J), models considered in [421 E3]> an d the MIT bag model itself. We neglect their 
divergent when T < 2~ 1//4 Td ec - Line (2) denotes the hard-wall model with the spherical 
horizon in Eq. (jSJ). Line (3) indicate the soft- wall model Eq. (JHj) shows. Line (4) 

indicate the ten-dimensional "AdS/QCD cousin" model in Eq. (J7J). Line (5) denotes the 
fuzzy bag model result for comparison with Line (4). Line (6) is for the Giirsoy et al. 
model given in [51]. Line (7) is calculated by the phenomenological model of |22j, with 
a scalar potential V((f)) = [—12 cosh (-\/7/120) + 2(j) 2 ]/L 2 . The thin gray lines are the 
p4-action result [53], in which the solid line indicates the pure glue case, the dashed line 
for the (2 + 1) flavor case, and the dashed-dotted line for the 3 flavor case. The points 
are calculated by the lattice methods with almost physical quark masses [54], where small 
solid bullets for N T = 4 case and solid squares for N T = 6 case. The small dark region 
near the critical temperature is enlarged and shown in the top-right corner, where the 
triangle-like shape formed by some line segments shows clearly the multi- valued nature of 
Line (6) and (7). 
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Figure 2: The entropy density s = —df/dT for various models compares to the free gas 
case. In fact, in the large- N c limit, we have s q oc iV c 2 and Sh oc N®, hence Sh = 0. The 
notations are as in Fig. [TJ The thick lines are models (1^4^5^7^2^3^6) 
respectively seeing from the arrow direction. 



26 



0.4 




-0.1 ^ '— ' ■— ' '— 1 ■— ' '— ' !— ' '— ' 

12 3 4 

T/r dec 

Figure 3: The square of sound speed = d log T/d logs for various models compares to 
the free gas case. For the ideal gas case, or the strongly coupled Af = 4 SYM theory- 
indicated in Eq. (JT]), we have Cg = 1/3. The notations are as in Fig. HJ The thick lines are 
models (1— >4-^5^7^6^2^3) respectively seeing from the arrow direction. 
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Figure 4: The normalization of the latent heat for pure glue fields of various models. In 
this case, we have (7t 2 /15)(As/ssb) = (Lh/T^ ec )/(N^ — 1). The left part is calculated for 
various models as tagged in Fig. HJ The right part shows that lattice results for iV c = 3, 4, 6, 
and 8 [37] . which suggest that the phase transition is second order for N c = 2, weakly first 
order for N c = 3, and robustly first order for N c > 4. The black error bars are for L t = 5, 
and the gray ones for L t = 6 and 8. The fitting line for L t = 5 case is informal, but the 
extend to iV c — > oo case can be guess in this fitting. Notice that the original MIT bag 
model has the value 7r 2 /15 ~ 0.658 in this figure. 
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Figure 5: The prefactor (K/27r)f2 in the nucleation rate formula. The thin gray dashed 
and dashed-dotted lines on the top are for dimensional values T c 4 and T 4 respectively. The 
thin gray solid line using the same parameters as in |lU4j . is shown for comparison reasons. 
Its value seems much larger than all other cases, mainly because it uses a rather large 
u = 50 MeV/fm 2 (although other parameters also affect the curve); however, we choose a 
rather small value of a = 0.02 T d 3 ec ~ 3.64 MeV/fm 2 for T dec = 192 MeV [105] in all other 
estimations. The gray solid bullet and square lines are for the pure gluon SU(3) lattice 
result L h = 1.4 T^ ec and o = 0.02 T| ec [37J. Nevertheless, when calculating the effectively 
massless degrees of freedom, we also count the fermionic contributions. The difference is 
the former case uses the perturbative result r/ q ~ 1.12 T 3 /a 2 log (l/o; s ) and a s ~ 0.23, but 
the latter case uses the AdS/CFT result r] q = s q /47r. The thick color lines are for models 
discussed above. Seeing from the arrow direction, they are models (7— > 2 — > 3 ^ 5 ^ 
4 — > 1) respectively. The process of scaling those large- N c theories to real QCD, and the 
rationality of that scaling, are discussed in the main text. The shear viscosity of models 
(1) and (4) are evaluated by [77]; while for all other cases, we choose r/ q = s q /47r. The 
bulk viscosities are chosen by the relation ( q /r] q = 2(1/3 — c 2 ) of Eq. ffl~5l) . and the shadow 
regions show the differences between them and the ( q = cases. ( q of model (7) can be 
calculated from more sophistical numerical results given by [221 121] if needed; however, we 
deal with it similarly with others for simplification. The black dotted line near the bottom 
is for the original MIT bag model with i] q = (s q — Sh)/47r. We choose the correlation length 
£ q = 0.48 (Tdec /T) fm [106j from lattice result for all our estimations, except the thin gray 
solid comparison line; in the gravity side, a lower limit of £ q is given by [67J. 
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Figure 6: The supercooling scale A = 1 — T f /T dec depends on the surface tension a for 
various models, which is estimated by Eq. (1231) . The notations are as in Fig. [5j except 
the shadow regions around the lines show the difference between Eq. (1231) and the rough 
criterion V ~ The two gray vertical dashed lines are marked for o = 0.02 T| ec and 

o = 0.2 Tj cc , which are chosen as typical parameters in Fig. [H]and[IUJ It can be seen that 
the supercooling scale A is really unsensitive to the method we estimate it, even in the 
small a regions where d nnc -C g?h- The thick lines are models (7^2^3^5^4— > 
1) respectively seeing from the arrow direction. Although t> s h can be calculated accurately 
by [90], we choose t> s h = c s for simplification, where the differences between them are 
imperceptible. 
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Figure 7: The mean nucleation distance d nuc depends on various surface tension a, esti- 
mated by Eq. ( 1241) . The notations are as in Fig. [5j The thick lines are models (7^2 
^3^5^4^1) respectively seeing from the arrow direction. Although the terminal 
point "*" marked for model (7) is factual, the terminal point "o" marked for model (2) is 
the numerical limit of our calculation. A maximum a exists for the maximum expected 
supercooling scale to be achieved; as Lh = for A = A< in model (2), d nuc — > oo when a 
tending towards this limit. 




Figure 8: A Hawking-Page phase transition [21] should always have a minimum tempera- 
ture T min , below which the high temperature phase cannot exist. This minimum tempera- 
ture is intrinsic, rather than caused by impurities or perturbations in the old supercooling 
scenarios. The long curved arrows show the behavior of the system from high temperature 
to low temperature phase, if no supercooling happens. 
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Figure 9: The mean nucleation distance d nvLC depends on the latent heat for the Gubser 
et al. model [22]. The gray vertical dashed line is = 2.67 T dec deduced from the potential 
V((f)) = [—12 cosh (a/7/120) + 2(f) 2 ]/ L 2 in the formal estimations. The three thick lines are 
for a = 0.2 T d 3 ec , 0.02 T d 3 ec and 0.002 T d 3 ec (from up down), respectively. 



32 



1.0 



1.5 2.0 

LiJTiicc 4 



3.0 



Figure 10: The various supercooling scales depend on the latent heat L h for the Gubser et 
al. model [22] . The black dotted curve from the top-right corner to the bottom- left corner 
is the maximum supercooling scale A<; hence the shadow region above it, is forbidden by 
that model. The actual supercooling scales A = 1 — T f /T dec calculated by Eq. (p^|) are 
denoted by the thick solid lines, which are for a = 0.2 T| ec , 0.02 T| ec and 0.002 T| ec (from 
up down) respectively. The dotted lines a little below them constrain the phase transition 
to be completed, which are roughly calculated by T ~ that is, d nvLC ~ d^. For some 

particular a, d nuc can easily be much larger, providing that the latent heat is small 
enough that the maximum supercooling A< is saturated. However, it is unlikely that the 
larger d nuc can help us understanding the formation of quark nuggets or the inhomogeneous 
initial conditions of the big-bang nucleosynthesis, because the parameter L h should be fine 
tuned. 
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